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Abstract 

We revisit the three-body problem in the framework of general relativity. The Newtonian N-body 
problem admits choreographic solutions, where a solution is called choreographic if every massive 
particles move periodically in a single closed orbit. One is a stable figure- eight orbit for a three- 
body system, which was found first by Moore (1993) and re-discovered with its existence proof by 
Chenciner and Montgomery (2000). In general relativity, however, the periastron shift prohibits 
a binary system from orbiting in a single closed curve. Therefore, it is unclear whether general 
relativistic effects admit a choreographic solution such as the figure eight. We carefully examine 
general relativistic corrections to initial conditions so that an orbit for a three-body system can be 
closed and a figure eight. This solution is still choreographic. This illustration suggests that the 
general relativistic N-body problem also may admit a certain class of choreographic solutions. 

PACS numbers: 95.10.Ce, 95.30.Sf, 45.50.Pk, 04.25.Nx 
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Introduction. — The three-body problem in the Newton gravity is one of classical problems 



in _ y „d P n yslcs (e,, I). „ 1765 , Euter found a collar sol» tion , and Lag^e 
found an equilateral triangle solution in 1772. It is impossible to describe all the solutions 
to the three-body problem even for the 1/r potential. In fact, Poincare proved that we 
cannot analytically obtain all the solutions, and the number of new solutions is increasing 
j^]. Therefore, the three-body problem remains unsettled even for the Newtonian gravity. 
The Newtonian N-body problem admits choreographic solutions, which attract increasing 
interests. Here, a solution is called choreographic in the celestial mechanics if every massive 
particles move periodically in a single closed orbit. In fact, a choreographic figure-eight 
solution to the three-body problem was found first by Moore [3] and re-discovered with its 
existence proof by Chenciner and Montgomery (4]. 

The theory of general relativity is currently the most successful gravitational theory 
describing the nature of space and time, and well confirmed by observations. Especially, it 
has passed "classical" tests, such as the deflection of light, the perihelion shift of Mercury 
and the Shapiro time delay, and also a systematic test using the remarkable binary pulsar 
"PSR 1913+16" |5|. It is worthwhile to examine the three-body (or more generally, N- 
body) problem in general relativity. N-body dynamics in the general relativistic gravity 
plays important roles in astrophysics. For instance, the formation of massive black holes in 
star clusters is tackled mostly by Newtonian N-body simulations (e.g, However, it is 
difficult to work out in general relativity compared with the Newtonian gravity, because the 
Einstein equation is much more complicated [3] (even for a two-body s yste m 8 L 9, 10, llj|). 



In addition, future space astrometric missions such as SIM and GAIA 12J, ll3|, llJ] require a 



general relativistic modeling of the solar system within the accuracy of a micro arc-second 



151 ] . Furthermore, a binary plus the third body were discussed also for perturbations of 
gravitational waves induced by the third body jltl 17 . Q, 19]. In this paper, we do not 



intend to solve the N-body problem in general relativity under a general situation. Instead, 
we shall focus on a choreographic solution. No choreographic solution has not been found to 
the general relativistic N-body problem so far. 

In a two-body system, the post-Newtonian corrections cause the periastron shift so that 
the binary system cannot orbit in a single closed curve |7j. As a result, it is unclear whether 
general relativistic perturbations admit a choreographic solution as the figure-eight. One 
may thus ask, "what happens for the figure-eight in the Einstein's gravity?" A specific 
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FIG. 1: A schematic figure for a binary orbit in general relativity. The orbit is not closed any 
more, because of the periastron shift. 



question may arise, such as "does the figure-eight cause periastron shift?", or "does the 
figure-eight make a transition to an open orbit in the general relativistic gravity?" The 
purpose of this paper is to answer these questions by carefully examining general relativistic 
effects to initial conditions for being a choreographic solution. 

This paper is organized as follows. First, we briefly summarize the choreographic figure- 
eight solution in the Newton gravity. Next, we analytically examine initial conditions and 
numerically solve the Einstein-Infeld- Hoffman equation of motion in order to obtain a choreo- 
graphic solution in general relativity. Throughout this paper, we take the units of G = c = 1. 
Newtonian choreographic solution. — As mentioned above, it is impossible to describe all 
the solutions to the three-body problem even for the 1/r potential. The simplest periodic 
solutions for this problem were discovered by Euler (1765) and by Lagrange (1772). The 
Euler's solution is a collinear solution, in which the masses are collinear at every instant with 
the same ratios of their distances. The Lagrange's one is an equilateral triangle solution in 
which each mass moves in an ellipse in such a way that the triangle formed by the three 
bodies revolves. Built out of Keplerian ellipses, they are the only explicit solutions. In 
these solutions, each mass moves on an ellipse. A choreographic solution for which three 
bodies move periodically in a single figure-eight orbit was found first by Moore by numerical 
computations 3J. The existence of such a figure-eight orbit was proven by Chenciner and 
Montgomery 4fl. This solution is stable in the Newtonian gravity 2fJ, |21[. The figure-eight 
seems unique up to scaling and rotation according to all numerical investigations, and at 
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the end its unicity has been recently proven 22j. Furthermore, it is shown numeri cally 



that fourth, sixth or eighth order polynomial cannot express the figure-eight solution [21]. 
Nevertheless, no analytic expression in closed forms for the figure-eight trajectory has been 
found up to now. Therefore, in this paper, we numerically prepare the figure-eight orbit. 

For simplicity, we assume a three-body system with each mass equal to m. Without loss 
of the generality, the orbital plane is taken as the x — y plane. The position of each mass 
(jtia) is denoted by (xa,Va) for A = 1,2,3. Figure [2] shows the figure-eight orbit, where we 
take the initial condition as £ = (xi,yi) = (— x 2 , — y 2 ) = (97.00,-24.31), (x 3 ,y 3 ) = (0,0) 
and V Newton = (£3,3/3) = {-2x u -2 yi ) = (-2x 2 ,-2y 2 ) = (-0.09324,-0.08647) , where 
a dot denotes the time derivative [2lJ. When one mass arrives at the knot (center) of the 
figure-eight, £ = \£\ is a half of the separation between the remaining two masses. It is 
convenient to use £ instead of a distance between the knot and the apoapsis, because the 
inertial moment is expressed simply as 2ml 2 . The orbital period is estimated as ^Newton = 
6.326m" 1 / 2 £ 3 / 2 « l^{M Q /mf/ 2 {£/R & f/' 2 sec, where M and R Q are the solar mass and 
radius, respectively. Obviously this system has no Killing vector as seen in Fig. El Here, we 
should note that £ is taken as 100, while it is the unity in the previous works. This is because 
we will treat the post-Newtonian correction in terms of the ratio between the mass and the 
separation such as £. In our case, the ratio m/£is 0.01, that is, the post-Newtonian correction 
becomes about one percent. In the equation of motion, the second post-Newtonian (2PN) 
corrections of the order of (m/£) 3 can be safely neglected, if m/l is very small, say, 10 -8 . In 
this case, however, the Newtonian and relativistic orbits will be indistinguishable. In order 
to demonstrate the difference between the two orbits, we choose m/£ as 0.01, for which the 
first post-Newtonian (1PN) terms are several dozens times larger than 2PN ones. In our 
computations, which are not long-time integrations over, say, thousands orbital periods, we 
can assume that 1PN terms are enough to bring major relativistic effects. 
Post- Newtonian figure-eight. — In the previous part, the motion of massive bodies follows the 
Newtonian equation of motion. In order to include the dominant part of general relativistic 
effects, we take account of the terms at the first post-Newtonian order. Namely, the motion 
of the massive bodies obeys the Einstein- Infeld-Hoffman (EIH) equation of motion [7J]. The 
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EIH equation is derived also from the first post-Newtonian Lagrangian as 23] 
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Figure [2] shows an orbit of a body starting at the Newtonian initial condition described 
above. In Fig. [2} a figure-eight orbit does not seem to survive at the 1PN order. However, 
this is not the case. We should note that the initial condition at the 1PN order does not 
necessarily coincide with that for the Newtonian gravity. We will thus carefully examine the 
initial condition by taking account of 1PN corrections. For this purpose, we assume that 
both the linear and angular momenta are zero (i.e. P = and L — 0). 

We should remember Vi = v 2 = —(v 3 )/2 for the Newtonian figure eight, for which both 
the total linear and angular momenta are zero. The changes in V\ and v 2 are expressed by 
using two vectors v% and £, which are linearly independent. 

Hence, the initial velocity of each mass is parameterized as 

777 

v ± = kV + t-{V-t)l, (5) 

777 

v * = kV + t-{V-t)l, (6) 
v 3 = V, (7) 

where k is expressed as 

k = - 1 - + a\Vf + ^. (8) 

Here, the 1PN terms have either \V\ 2 or m/t. If P = and L = 0, there is no need of | V\ 2 
in front of I in Eqs. (JSj) and (jSJ) as shown below. 
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FIG. 2: Figure-eights starting at the Newtonian initial condition. The solid curve denotes a 
figure-eight orbit in the Newtonian gravity. The dashed curve denotes a trajectory of one mass 
following the EIH equation of motion under the same Newtonian initial condition. 



The linear and angular momenta are calculated from the first post-Newtonian Lagrangian 



23l |. Here, we impose the condition of P = and L = at 1PN order. Then, we determine 



the 1PN coefficients as 



a = — - 



16' 



(9) 
(10) 

(11) 



Ut to this point, V is arbitrary. Next, we determine V . 

The initial velocity of the particles can be different from that for the Newtonian gravity. 
The post-Newtonian effects affect both the magnitude and direction of the velocity. There- 
fore, by using two linearly independent vectors, £ and V Newton, we parameterize the initial 
velocity as 

/ 777 \ 777 £ ( £\ 

V = (1 + 5 — J V N ewton + V—j ( V New ton ' J J , (12) 

where it is sufficient to express 1PN corrections in terms of either m/£ or |VN C wton| 2 in 
numerical computations, though both are necessary for analytic calculations of P = and 



L = 0. For convenience's sake, we choose m/£ in Eq. (TP2]) . 

By numerically performing trial and error iterations until achieving a periodic orbit, we 
find out 



n 



s 



-3.3, 



-3.7. 



(13) 
(14) 



Here, the iterative computations are done until we find the values of 5 and i] for which the 
three masses simultaneously return to their initial positions. Our procedure is as follows. If 
and only if one particle returns to the neighborhood of its initial position (i.e. the origin 
for the particle labeled by 3) within the positional deviation of 0.01, we measure how far 
the remaining two particles are from their initial positions at the same moment when the 
particle is closest to the initial position. For Eqs. (TT3"j) and (TH|) . the sum of the square 
distances is minimized as approximately 0.1. This is sufficient for i = 100 and mjl = 0.01, 
because expected positional shifts after one cycle are of the order of the unity or more. For 
instance, such a shift exceeds ten for a = f3 = !;=5 = r] = 0m Fig. [2j 

We integrate the motion over ten cycles to confirm the periodicity. After ten periods, 
the found solution comes to the same point within the deviation of ±1. The numerical 
computation gives the orbital period as 



The relativistic figure eight appears to be stable, in the sense that we recognize 'eight-like' 
orbits that resemble figure eight for slightly different values of S and r]. It is a future subject 
to analyze the long-time stability of the relativistic figure eight. 

Figure [3] shows that a figure-eight orbit is still closed even after including the dominant 
general relativistic effects. In Fig. [3l we can recognize an asymmetric difference between the 
Newtonian figure-eight orbit and the general relativistic (GR) one. The deviation is partly 
fiducial, because the principal axes of the GR figure-eight orbit are not along the x and y 
axes. That is, £, which defines the direction of the initial position of a particle with respect to 
the principal axes, changes slightly at the 1PN order. The axes are inclined by 0.012 radian 
with respect to those of the Newtonian figure-eight orbit. Now, we choose the X~clXlS clS 
the principal axis for both the Newtonian figure-eight orbit and the general relativistic one. 
After choosing the principal axes, Fig. H] shows a general relativistic choreographic solution 




(15) 
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FIG. 3: Figure-eight orbits. The solid curve denotes a figure-eight orbit in the Newtonian gravity. 
The dashed curve denotes a figure-eight orbit at the 1PN order of general relativity. 

at the first post-Newtonian order. The solution recovers line symmetry with respect to the 
x and y axes. There are no significant differences in the velocity between the Newtonian 
and GR figure-eight orbits. One may notice that the 2PN terms are neglected. It would be 
safer to choose mjl for instance as 10 -8 and then on the figures to exaggerate the differences 
between the Newtonian and post-Newtonian solutions in order to make them visible. 

Finally, we mention the possibility of three-body systems in a choreographic orbit such 
as a figure-eight. As a new outcome of binary-binary scattering, the figure-eight orbit was 
discussed for presenting a way of detecting such an orbit in numerical computations 24]. 
According to the numerical result, the probability of the formation of figure-eight orbits is 
a tiny fraction of one percent. The gravitational waves emitted by the figure-eight have 



been recently studied by assuming the motion in the Newton gravity [25|. By evaluating the 
radiation reaction time scale, it is shown also that figure-eight sources emitting gravitational 
waves may be too rare to detect. 

Conclusion. — We obtained a general relativistic initial condition for being a figure-eight 
orbit. This condition provides the first choreographic solution taking account of the post- 
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FIG. 4: Relativistic figure eight. The principal axes of the orbit are chosen as x and y axes. 

Newtonian corrections. It is interesting to include higher post-Newtonian corrections, es- 
pecially 2.5PN effects in order to elucidate the backreaction on the evolution of the orbit 
due to the gravitational waves emission at the 2.5PN order. If the system is secularly stable 
against the gravitational radiation, one might see probably a shrinking {£ < 0) figure-eight 
orbit as a consequence of a decrease in the total energy (E < 0). This speculation will 
be confirmed or rejected in future. It may be important also to look for other relativistic 
choreographic solutions for a system including four or more masses. It is possible that some 
of Newtonian choreographic solutions are prohibited by general relativistic effects. Further 
investigations along these lines will allow us to probe many-body dynamics in the Einstein 
gravity. 
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